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Abstract 

Hamiltonian systems with linearly dependent constraints (irregular systems), are classi- 
. tied according to their behavior in the vicinity of the constraint surface. For these systems, 

the standard Dirac procedure is not directly applicable. However, Dirac's treatment can be 
slightly modified to obtain, in some cases, a Hamiltonian description completely equivalent 
to the Lagrangian one. A recipe to deal with the different cases is provided, along with a 
CO , few pedagogical examples. 
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1 Introduction 



Dirac's Hamiltonian analysis provides a systematic method for finding the gauge symmetries 
present in a theory. The analysis identifies and classifies the constraints, which are local func- 
tions of the phase space coordinates. Consistency requires that the constraints be preserved in 
during the evolution (for the review of the Hamiltonian analysis see Ref.pQ-(l]). However, if the 
constraints are not functionally independent, then Dirac's procedure is not applicable. The test 
of functional independence are the so-called regularity conditions, and those systems which fail 
the test are said to be irregular. 

Irregular systems are not necessarily intractable or exotic. A simple example is a relativistic 
massless particle (p^p^ = 0), which is irregular at the origin of momentum space (p^ = 0). This 
point in phase space is exceptional, as it is unclear whether this would be an observable state 
for a photon, say. On the other hand, we know the configuration p^ = to be a very important 
one: the ground state. There are other physical circumstances in which regularity is violated, 
and not only for isolated states but on large portions of the region in phase space where the 
system evolves. Chern-Simons theories in 2n + 1 spacetime dimensions are examples where, for 
some initial configurations, regularity can fail at all times and one is forced to live with this 
problem. 

Here we discuss the possible ways in which the constraints can fail the test of functional 
independence, and how the Hamiltonian treatment of Dirac must be modified in each case. 
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2 Regularity Conditions 



If we call z % = (q,p) (i = l,...,2n) the coordinates in the phase space T, the constraints 
<p r {z) ~ (r = 1, ...R) define the constraint surface £ given by 



£ = {z G r | r (z) =0 (r = 1, . . . , R) (R < 2n)} 



(1) 



Regularity Conditions (RCs): The constraints (j) r ~ are regular if and only if their 
small variations 5(jf , evaluated on S, are R linearly independent functions of bz l . 



To first order in 5z l , the variation of the constraints are 5(j) r 
Consequently, the RCs can also be defined as [2]: 

The set of constraints 4> r ~ is regular if and only if the Jacobian J\ = 
rank: K(J) = R. 



JJ6z\ where JT = 



has maximal 



A simple classical mechanical example of functionally dependent constraints occurs in a 2- 
dimensional phase space with coordinates (q,p) and constraints (j) 1 = q ~ and (p 2 = pq ~ 0. 



In this case, J 



1 p 
q 



and 5? 



d(q,p) 



q=0 



fail the test of regularity. Consider the constraint 

" 2? " 



In this case, J 







1. A system of just one constraint can also 
b = q 2 ~ in a 2-dimensional phase space. 
and 3ft (J) = 0. The same problem occurs with the constraint 



q 2 =0 



q 7 ~ 0, which has a zero of seventh order at the constraint surface, or with any other constraint 



3 Classification of Irregular Constraints 

Irregular constraints can be classified according to their approximate behavior near the surface 
S. 

A. Linear constraints. The Jacobian has constant, non-maximal rank throughout E, and 

4> r = JT (V - A w 0, 3ft( J) =R' <R. (2) 

These are regular systems in disguise. Regularity fails simply because R — R' constraints are 
redundant and should be discarded. The regular system gives the correct description. 

B. Multilinear constraints. In the vicinity of S, the constraints are of the form 

= ^5 n ... im (/ 1 -^)---(^-z^) f»0 (m<2n), (3) 

where the coefficients Si u „i m vanish if any two indices are equal. Thus, <f> has simple zeros on 
surfaces of dimension 2n — 1, and zeros of higher order occur at the intersections of these surfaces. 
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The RCs fail at the points of intersection, where 4> has multiple zeros. At those intersections, (ft 
can be replaced by the equivalent 1 set of constraints, 

ip 1 = z h - z h w 0, • • • , <p m = z lm - z lm » 0. (4) 

At the intersections, the set {if 1 , f m } is regular and therefore this provides a recipe for 
substituting the irregular multilinear constraint <j) by a regular set of linear constraints. For 
example, the constraint 4> = qp « is irregular at (0, 0) because it admits a linear approximation 
everywhere except at this point. Replacement <j> by the linear constraints {q ss 0,p ~ 0} at (0, 0) 
regularizes the system. 

C. Higher order constraints. In the vicinity of S, the constraints do not possess a linear 
approximation: 

= C(z s -z s ) fc «O (jfe>l). (5) 

The Jacobian vanishes everywhere on the constraint surface. Naively, it would seem possible to 
choose z s — z s ~ as an equivalent regular constraint, but it turns out that this could change 
the dynamics of original theory, as we show below. 

These three classes are generic and, in general, there can be combinations of these three 
types occurring simultaneously near a constraint surface. By selecting a sufficiently small neigh- 
borhood of a point on S, one can always expect to be in one and only one of the three situations 
just described. 

In correspondence with the three generic cases in which regularity can fail, the nature of the 
constraint surface S falls into one of three categories: 

A. The RCs are satisfied on the whole constraint surface. These are regular systems, either 
desguised or not. 

B. The RCs fail on a submanifold of S: 5? (J) = R except on a submanifold So C S where 
sft(J) = R> < R . 

C. The RCs fail everywhere on S: J has constant lower than R rank on S. 

The first case is treated in the standard texts and will not be further discussed here. 

In the second case, the constraint surface can be decomposed into two non-empty submani- 
folds, So and S# such that S = So U S^ and So H S^ is empty. Then, the rank of the Jacobian 
jumps from 3ft (J) = R on Sr, to 3? (J) = R' on So- As mentioned above, in this case it is 
possible to replace (f> at So by a set of regular constraints {(p 1 « 0, • • • , (p m « 0} which regularize 
the system at So- 

The important question is how to proceed in the third case. If the RCs fail everywhere on 
S, the previous approach is not applicable, because there is no guarantee that the resulting 
Hamiltonian dynamics will be equivalent to that of the original Lagrangian system. 

This can be seen in the example of Lagrangian in a three-dimensional configuration space, 

L (x,y, z) = xz + yz 2 . (6) 
x Two sets of constraints are said to be equivalent if they define the same constraint surface S. 
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The general solution of Euler-Lagrange equations describes a system with one degree of freedom 
- a free particle, whose time evolution x = po t + xq is determined by two initial conditions po 
and xq. The remaining fields are y(t), a Lagrange multiplier with indeterminate evolution, and 
z(t), with trivial evolution, z(t) = 0. 

The Hamiltonian approach gives just two (first class) constraints, p y 0, z 2 (R = 2), 
and one degree of freedom, as expected from the Lagrangian approach. However, this is only 
superficially correct, because the Jacobian has rank 1, and not 2. This is because the constraint 
function (ft = z 2 has no linear approximation at z = 0. 

On the other hand, if we naively take z ~ 0, which is regular and equivalent to z 2 ~ 0, 
then the Hamiltonian analysis generates three first class constraints p y ~ 0, z ~ 0, p x ~ 0, 
which leave zero physical degrees of freedom. This result is not consistent with the Lagrangian 
description where there is one degree of freedom. 

As seen in the previous example, even if two constraints are equivalent in defining the same 
constraint surface, they may yield different dynamics and should be treated more carefully. 

Suppose (ft ~ is a constraint of the form @, which is equivalent to the regular constraint 

X = 1/fc « 0. (7) 

The question is whether x ~ gives also the correct dynamics. The answer to this question 
depends on whether x is a first °r second class constraint. Namely, it makes a difference whether 
the linear constraint x can generate a transformation in phase space that leaves the Hamiltonian 
action unchanged or not . As shown in jS], if the linearized constraint x is second class, then 
it is not only geometrically equivalent to (ft in the sense that it defines the same S, but the 
substitution also yields the same dynamical description as the Lagrangian approach. On the 
other hand, if x is first class, then the subtitution generates a system whose dynamics is different 
from the one obtained from the Euler-Lagrange equations. 

4 Conclusions 

The recipe for treating the non-regular constraints is: 

• Every linear or multi-linear set of constraints can be exchanged by an equivalent regular 
set. It allows to carry out Dirac's procedure in the standard way. 

• A higher order constraint (ft = C {z s — z s ) k ~ 0, can be exchanged by the equivalent linear 
constraint x = ~ 0- If X is a second class constraint, the dynamics of the new system 
is equivalent to the Lagrangian one. If x is first class, the substitution yields a system 
which is not dynamically equivalent to the Lagrangian one. In this latter case, one should 
view the original Lagrangian as an incomplete, if not a totally inconsistent description for 
a dynamical system. 
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